Local electronic structure near oxygen dopants in BSCCO-2212: 
a window on the high-T c pair mechanism? 
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The cuprate material Bi2Sr2CaCu208 (BSCCO-2212) is believed to be doped by a combination 
of cation switching and excess oxygen. The interstitial oxygen dopants are of particular interest 
because scanning tunnelling microscopy (STM) experiments have shown that they are positively 
correlated with the local value of the superconducting gap, and calculations suggest that the fun- 
damental attraction between electrons is modulated locally. In this work, we use density functional 
theory to try to ascertain which locations in the crystal are energetically most favorable for the O 
dopant atoms, and how the surrounding cage of atoms deforms. Our results provide support for the 
identification of STM resonances at -leV with dopant interstitial O atoms, and show how the local 
electronic structure is modified nearby. 

PACS numbers: 74.25.Bt,74.25.Jb,74.40.+k 



STM measurements of impurities and other inhomo- 
geneities in the cuprates have opened a new window on 
high temperature superconductivity and raised a host of 
new questions about the way these impurities interact 
with their environment 0, 0, Q • It has been traditionally 
assumed that an impurity acts as a localized screened 
Coulomb potential, and that its principal effect is the 
modification of quasiparticle wavefunctions nearby. This 
effect is large and observable as a consequence of the d- 
wave symmetry of the superconducting stateQ. Within 
such models, the superconducting order parameter is sup- 
pressed around the impurity, but this phenomenon is not 
usually essential for qualitative predictions. Recently, a 
phenomenological analysis Q of STM experiments imag- 
ing interstitial oxygen atoms 6] suggested that such im- 
purities might have a much more striking effect, to wit, 
the local modulation of the electronic pair interaction, 
leading to large amplitude modulations of the supercon- 
ducting order parameter. If this is true, an understand- 
ing of the local changes in electronic structure and cou- 
plings to collective excitations of the material might tell 
us which aspects are of crucial importance to the pairing 
interaction itself. 

Naturally occurring impurities which dope the Cu02 
planes in as-grown P^S^CaC^Os crystals 7j include 
both the roughly 3% excess Bi substituting on the Sr 
sites, as well as oxygen interstitials, whose concentration 
depends on the annealing sequence and determines the 
net doping of the sample. Both sets of dopant atoms 
have recently been imaged by STM@, 0. In the case 
of the O interstitials, which appear as bright spots on 
the Bi2Sr2CaCu208 surface at a bias of -960mV, a re- 
markable set of correlations was established by McElroy 
et al.[|| between the positions of the impurities and the 
magnitude of the local superconducting gap (as defined 
by the position of the coherence peak in the conductance 
signal), as well as the magnitude of the LDOS at differ- 
ent energies. Local doping disorder was therefore claimed 



to be responsible for the nanoscale inhomogeneities ob- 
served earlier in the same material \wL IllL Il3| . but 
the positive dopant-gap correlation observed contradicted 
earlier theoretical proposals along these lines [Til ITU Il6 | . 
Nunner et al.Q then proposed that the oxygen dopant 
increased the pairing strength g between electrons locally, 
and a simple model based on this ansatz was shown to 
explain all experimentally reported correlations at op- 
timal doping. The results suggested that the range of 
the modulation was of order a Cu-Cu distance. While 
such atomic-scale enhancements of pairing are outside 
the framework of standard BCS-type theoretical treat- 
ments of impurities in superconductors, they have occa- 
sionally been discussed in earlier treatments of the ef- 
fects of impurities and twin boundaries on conventional 
superconductors |17| . 

If the theoretical proposal is correct, understanding 
where the O interstitials sit and how they distort the 
cage of atoms around them could be an important clue to 
the ingredients of the electronic pair mechanism in these 
materials. Even if not, the details of the local chemistry 
around interstitials can be important to understand the 
doping process in the cuprates. At present, there is no 
consensus as to the exact position of the defects imaged 
by STM, their chemical identity, nor the nature of the 
states they induce at -960 meV into which one tunnels 
in the experiment. We have therefore undertaken den- 
sity functional calculations of BSCCO-2212 including O 
atoms placed in the most natural interstitial locations be- 
tween the BiO and SrO, or between the SrO and Cu02 
planes. We regard locations between other layers as un- 
likely, since they are too far from the BiO layer, where the 
BSCCO-2212 crystal typically cleaves, to be the source 
of the strong STM signal observed. They are allowed to 
relax to lower the overall energy of the system, and ul- 
timately several possible stable and metastable positions 
are found. Only one of these, close to but not identical 
to the position found by statistical analysis of the STM 
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1 


1.09,1.09,-0.93 


-6.56 


2 


2.11,0.00,-0.38 


-3.91 


3 


2.60,0.00, -0.38 


4.38 


4 


2.60,0.00,-4.81 


14.97 



TABLE I: Stable/metastable positions of interstitial found. 
Coordinates are relative to central Bi atom in Fig. Energies 
are per O interstitial measured relative to crystal without 
interstitial. 

data, represents a clear minimum of the energy. We an- 
alyze the structural and electronic changes thereby pro- 
duced in some detail. 

The most important results include a displacement 
of the nearby apical (stoichiometric) oxygens from their 
naive positions on the Bi-Cu axis and the creation of 
states centered primarily around -leV, not only associ- 
ated with the interstitial O itself, but with the apical 
oxygens and nearby Cu sites. It is the unhybridized 
dopant 02p z states which give rise to the resonance seen 
by STM. It is remarkable that an ab initio calculation in 
an extremely complicated material is able to so quanti- 
tatively reproduce the energetics of such states, and we 
believe the prospects are good, in light of this success, 
for wider applicability of the technique to study other 
impurity problems in the cuprates. 

Method. In this work, DFT calculations have been 
used to determine all structural, energetic and elec- 
tronic results. The Kohn-Sham equations are solved 
self-consistently in a plane- wave basis set, in conjunction 
with Vanderbilt ultrasoft pseudopotentialsjl^lljj, which 
describe the electron-ion interaction, as implemented in 
the Vienna ab initio simulation program (VASP)[2(|. 
Exchange and correlation are described by LDA. We 
use the exchange-correlation functional determined by 
Ceperly and Alder |2*lj and parameterized by Perdew and 
Zunger 22] . Details of the application of the method to 
impurity problems in Bi2Sr2CaCu2 0s have been given in 
|23j . It is known that LDA fails to give the correct band 
structure for cuprates in general, and in particular fails 
to describe the metallic/insulating character of the par- 
ent compounds. On the other hand, the bands and Fermi 
surface given by LDA agree relatively well with ARPES 
results on optimally doped BSCCO-2212, with minor ex- 
ceptions (see below). Furthermore, as argued in [23], we 
are primarily interested in high-energy, localized impu- 
rity states which should be little affected by the strong 
electronic correlations which renormalize the band states 
near the Fermi surface. 

The BSCCO-2212 surface is nearly always exposed 
at the cleaved BiO layer. For the calculation of the 
BSCCO surface in the presence of impurities, we there- 
fore use surface unit cells (super cells) with dimension of 
(2y2a x 2\/2a)R45 that contain 8 primitive unit cells in 
the xy plane. There is one O interstitial added to the sto- 
ichiometric sample of 120 atoms. The vacuum between 




FIG. 1: (Color online) Atomic displacements due to the pres- 
ence of O interstitial at position 1 in table [I] Atomic species 
are color coded: Bi (blue), stoichiometric O (yellow), Cu 
(green), Sr(light blue), Ca(aquamarine), O interstitial (red). 
Arrows indicate rough measure of displacement of atoms from 
their positions without interstitial. 



neighboring slabs is about 15 A thick. 

Structural relaxation. We compute the structure and 
electronic wavefunctions for the BSCCO-2212 surface 
first without and then with an O interstitial. In the lat- 
ter case, the O is positioned somewhere in the lattice and 
the entire system is allowed to relax to minimize the total 
energy. Several metastable positions were found, summa- 
rized in Table [i] It is seen that only one minimum thus 
found is a plausible location (position 1 in the table) for 
such a defect; others are higher in energy by at least 2.5 
cV. Position 1 is located in the plane between 2 Bi atoms 
in the BiO layer and 2 stoichiometric oxygens in the SrO 
layer below, as shown in Fig. ^ This position is dis- 
placed by about 0.5A and has a different symmetry from 
that determined experimentally by the STM group[f| by 
doing a statistical analysis of their -960meV signal. This 
may be attributable to a combination of dopant cluster- 
ing effects and the distortions due to the incommensurate 
supermodulation in this material. 

The displacements of atoms near the O impurity in 
position 1 are also shown schematically in Figure ^ and 
summarized in Table [H] The largest displacements are 
in the BiO and SrO planes. In particular, we note that 
the so-called apical oxygen in the SrO plane is moved by 
nearly half an A away from the Bi-Cu axis, leading to 
an effective tilting of the semi-octahedron of stoichiomet- 
ric oxygens away from this axis. Crudely speaking, this 
tilting is of the LTT (low-T tetragonal) type, which has 
been associated with the tendency towards stripe forma- 
tion in the underdoped cuprates[24|. Finally, the Bi and 
O(Bi) atoms nearest the impurity in the BiO layer are 
also strongly pushed away by displacements of the same 
order. 

Changes in electronic structure. Within the same 
framework, we now calculate the electronic structure and 
changes thereto caused by the O interstitials. To com- 
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Atom 


r*(A) 


rf-n (A) 


|r/ - r, 


Bi 


0,0,0 


-0.34,-0.34,-0.14 


0.50 


Bi' 


2.59,2.59,0 


-0.17,-0.17,-0.08 


0.25 


O(Bi) 


2.59,0,0.05 


0.68,-0.31,0.02 


0.75 


Sr 


2.59,0,-2.76 


0.0,0.0,0.10 


0.10 


O(Sr) 


0,0,-2.05 


-0.41,-0.41,-0.12 


0.59 


O(Sr)' 


2.59,2.59,-2.05 


0.16,0.16,-0.02 


0.23 


Cu 


0,0,-4.60 


-0.02,-0.02,0.14 


0.14 



(A) 



TABLE II: Initial positions r$, displacements Tf—Ti, and dis- 
placement magnitude of various nearby atoms in the presence 
of interstitial O at position 1 (see Table 1). Atoms are those 
of given species closest to O interstitial, unless indicated by a 
' (second closest). 

pare directly to the STM measurements, we should cal- 
culate the local density of states (LDOS) a tip height 
distance (~ 3. 5-4 A) above the BiO plane. At present our 
numerical resolution is insufficient to address this ques- 
tion directly, but we can ask in which energy ranges new 
states are found, and roughly where in space they arc lo- 
calized. The partial density of states (PDOS) projected 
onto atomic species \x is defined as 

M e ) = $3«fc$3*(e-ei,k)|fa„,kl**,k>| 2 > (1) 

k i 

where Ml k are linear combinations of atomic wavefunc- 
tions consistent with the crystalline symmetry, and the 
^i,k are Bloch wavefunctions for the ith. band. Here a?k 
is the weight of each k point and e^k is the dispersion of 
the Bloch state. 

In Figure [3 we plot the PDOS for relevant atomic or- 
bitals in the absence of the O interstitial, together with 
the same quantities with the interstitial present. The 
spectra have been plotted relative to the absolute Fermi 
levels of 0.6527 eV and 0.4832 eV for the case without and 
with O interstitials, respectively. The PDOS and DOS 
for the crystal without interstitials are consistent with 
earlier LDA calculations for this material [25J. Upon ad- 
dition of the O interstitial, one obvious change is the shift 
of the Cu and O(Cu) spectra towards the Fermi level, 
consistent with the hole doping of the Cu02 plane by 
the interstitials. Another interesting feature is the shift 
of PDOS Fermi level weight on the Bi sites to above the 
Fermi level upon doping with O. All previous LDA- type 
calculations of stoichiometric Bi-2212 have found BiO- 
derived bands crossing the Fermi level, leading to BiO 
Fermi surface "pockets" |23. These features have never 
been observed by ARPES|27|. suggesting that the dop- 
ing by excess O may be necessary to obtain the correct 
Fermi surface pjjj. 

The primary effect of the interstitial O appears to be to 
contribute states in a range around -1 to -3.5 eV, with a 
sharp concentration of states around -1.0 eV itself. Since 
there are originally no Bi, O(Bi), or O(Sr) states in this 
range, the extra O apparently causes a significant change 
to the PDOS of these species, inducing peaks at -1.0 eV in 



all three cases. O interstitial states in the rest of the en- 
ergy regions do not appear to induce qualitative changes 
in the PDOS of atomic species near the surface, and may 
therefore be expected to not cause significant qualitative 
changes to the LDOS measured by STM. The only ex- 
pected strong signal of the presence of the O defects is 
therefore a concentration of (assumed) tunnelling states 
at -1.0 eV, nearly precisely where such a signal is found 
by STM 6] . This is a spectacular confirmation of the abil- 
ity of DFT to predict electronic structure near impurities 
with high accuracy in these systems. We note that the 
rough magnitude of the binding energy of these states can 
be understood by recognizing the similarity to the stoi- 
chiometric 02p derived states (found at ~ -1.5eV below 
the top of the Zhang-Rice band) in photoemission ex- 
periments on Sr2Cu02Cl2, which do not hybridize with 
Cu 3d states[2Jj. The peaks around -leV seen in the 
O(Cu) and Cu spectra in Figure 2 which do not change 
significantly with the addition of O dopants are of this 
type. 

We note finally the existence of a broad cluster of 
O(dopant) states centered around -1.7 eV which also ap- 
pear to induce similar features in O(Bi) and Bi, and could 
lead to significant LDOS change at the BSCCO surface. 
This is an independent prediction which can be checked 
by STM. 
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FIG. 2: (Color online) Projected densities of states from Eq. 
0for different atomic species before (black solid line) and after 
(red dashed line) addition of O interstitial. 

Local density of states near -leV. The PDOS is calcu- 
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lated for all atoms of a particular species and symmetry, 
but spatially resolved LDOS information is necessary to 
specify where these changes are taking place relative to 
the interstitial site. To this end, we focus on the states 
near -1 eV, and plot a contour of the LDOS 

p(e,r) = ^ Wk ^^(6-6 4 , k )|^ k (r)| 2 (2) 

k i 

integrated over a range of energies from -l.leV to -0.9eV 
in Figure |3J Inspection of various contours of this type 
leads to the conclusion that the states at -leV are dom- 
inated by unhybridized Cu02 plane states [23, and in- 
duced states localized near the impurity. In particular, 




FIG. 3: (Color online) Spatial distribution of LDOS integrat- 
ing over states near -leV. The surface of the sample is at the 
top of the figure. Atomic species are color coded: Bi (blue), 
O(yellow), Cu (green), Sr(light blue), Ca(aquamarine), O in- 
terstitial (red). 

the slightly misoriented 02p z state on the O interstitial is 
clearly visible. There is no other state in the vicinity close 
to the surface positioned and oriented to accomodate an 
electron transfer from an STM tip near the surface. We 
note the absence of any states localized on the Bi atoms, 
in apparent contradiction to the results shown in Fig. [21 



This is because the Bi radius is 2 times larger, such that 
the PDOS weight appears at smaller LDOS values than 
depicted here. 

The next fascinating change due to the O interstitial 
is the creation of new states on the closest Cu site which 
appear to be of 3d 3z 2_ 1 symmetry. This process is clearly 
mediated by the apical O(Sr), which also shows induced 
%Px,y symmetry states in this range. There is no obvi- 
ous hybridization between these states; rather, these ap- 
pear to be changes due to the displacements of the atoms 
in question, which transfer spectral weight from higher 
binding energy into this range, as seen in Fig. 

Conclusions. Using density functional theory meth- 
ods, we have shown that the excess O dopant atoms in 
the as grown BSCCO-2212 material are located with high 
probability between the BiO and SrO layers in the plane 
containing Bi and O(Sr) atoms. These interstitials are re- 
sponsible for the creation of a set of states narrowly local- 
ized near -leV, which we have investigated more closely. 
One unhybridized 2p state appears on the dopant itself, 
and is oriented in the z direction, strongly suggesting 
that this is precisely the state imaged in recent STM 
experiments ,6] at -960 meV. This result is a strong con- 
firmation of the ability of DFT to calculate high-energy 
impurity features even in strongly correlated systems, as 
emphasized in Ref. |23l |. In addition, new states in this 
range are induced on the strongly perturbed nearby api- 
cal oxygen site, and even on the closest Cu atom. This 
suggests that the pairing of two electrons in the Cu02 
plane in this high-temperature superconductor may in- 
deed be influenced by the presence of the nearby O inter- 
stitial, as suggested in Ref. [f|, and that information on 
the modulation of the pair interaction may be gleaned by 
"downfolding" the DFT information to determine how 
the local change in electronic structure correlates with 
the superconducting gap. 
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